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We study the thermodynamic properties of a small superconducting metallic grain using a quantum Monte 
Carlo method. The grain is described by the universal Hamiltonian, containing pairing and ferromagnetic ex- 
change correlations. In particular, we study how the thermodynamic signatures of pairing correlations are 
affected by the spin exchange interaction. We find the exchange interaction effects to be qualitatively different 
in the BCS and fluctuation-dominated regimes of pairing correlations. 
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I. INTRODUCTION 

The properties of conventional bulk superconductors are 
well described by the Bardeen-Cooper-Schrieffer(BCS) mean 
field theory. 1 The BCS theory is valid in the limit when the 
pairing gap A is much larger than the single-particle mean- 
level spacing 8. However, in small metallic grains, the dis- 
creteness of the spectrum is important and the mean-level 
spacing can be comparable or larger than the pairing gap. This 
is the fluctuation-dominated regime, in which BCS theory is 
no longer a good approximation. 

The reduced BCS Hamiltonian was used extensively to 
study the properties of small metallic grains. 2 "— It was found 
that pairing correlations in the crossover between the bulk 
BCS limit and the fluctuation-dominated regime manifest 
through the number-parity dependence of thermodynamic 
quantities such as the spin susceptibility 9 " 12 and the heat ca- 

111-3 

pacity»ii— 

However, the effective low-energy interaction between 
electrons in a metallic grain contains additional terms be- 
yond the reduced BCS Hamiltonian. Such residual interac- 
tions could have significant effects on the signatures of pairing 
correlations in a finite-size grain. Finding this effective inter- 
action is, in general, a difficult task. However, a remarkably 
simple effective Hamiltonian emerges in grains whose single- 
particle dynamics are chaotic or weakly diffusive (in the pres- 
ence of disorder) in the limit of a large Thouless conductance 
gr i 14 ' 15 In such grains the single-particle Hamiltonian of ~ gj 
levels around the Fermi energy is described by random ma- 
trix theory. 16,17 The randomness of the single -particle wave 
functions induce randomness into the corresponding electron- 
electron interaction matrix elements. These matrix elements 
can then be decomposed into their average and fluctuating 
parts. The average interaction is determined by symmetry 
considerations ; 14 ' 1 8 and includes, in addition to the classical 
charging energy, a Cooper-channel BCS-like interaction and 
an exchange interaction that is proportional to the square of 
the total spin S of the grain. This average interaction together 



with the one-body Hamiltonian describe the so-called univer- 
sal Hamiltonian. 1415 Residual interaction terms are of the or- 
der 1/gr and can be ignored in the limit of large gr- 

Much work has been devoted to the understanding of pair- 
ing correlations in finite-size systems and, in particular, in 
small metallic grains. 8 Exchange correlations were also stud- 
ied extensively in semiconductor quantum dots, 19 " 22 where 
the pairing interaction is repulsive and can thus be ignored. 
Much less is known about the properties of a superconducting 
grain in the presence of both pairing and exchange correla- 
tions. 

In Ref. |23] we studied the phase diagram of the ground- 
state spin of a metallic grain that is described by the univer- 
sal Hamiltonian. The competition between superconductiv- 
ity and ferromagnetism leads to a narrow coexistence regime 
in the J s /h — A/8 plane. This regime can be broadened and 
tuned by an external Zeeman field. Signatures of this coex- 
istence were identified in the mesoscopic fluctuations of the 
conductance peak spacings and conductance peak heights in 
a metallic grain that is weakly coupled to leads. 24 Here we 
study the competition between pairing and exchange correla- 
tions in thermodynamic properties of the grain. In particular, 
we determine how the signatures of pairing correlations are af- 
fected by the spin-exchange interaction. Our studies cover the 
crossover from the fluctuation-dominated regime to the BCS 
regime. They are based on a quantum Monte Carlo method 
that is particularly suitable for the universal Hamiltonian. 

The outline of this paper is as follows: the model we use to 
describe the metallic grain (i.e., the universal Hamiltonian) is 
discussed in Sec. [Til while the quantum Monte Carlo method 
and in particular its application to the universal Hamiltonian is 
explained in Sec. [HI] Various thermodynamic properties are 
calculated in Sec. [TV] In particular, we discuss the thermal 
spin distributions (Sec. IIVA1) . the number of S = electron 
pairs (Sec. lIV Bl . the canonical pair gap (Sec lIVCb . the heat 
capacity (Sec. |IV Dl > and the spin susceptibility (Sec. |IV Ej . 
Our conclusions are given in Sec.lVl 
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H. THE MODEL 

The universal Hamiltonian of a metallic grain is given 



H 
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where d\„ are creation operators of electrons in spin- 
degenerate (a = ±) single-particle states with energy 
£/b N — Y^ka^la^kc * s tne particle-number operator, S = 
jUkaa' ^l c a c,a'C ka i is the total spin operator of the grain (a 

are Pauli matrices), and = Y^k^'k - 15 tne P a * r creat i° n 
operator in time-reversed (spin up/spin down) orbitals. Eq 
is the charging energy of the grain, while the parameters G 
and J s are the coupling constants in the Cooper channel and 
in the exchange channel, respectively. The universal Hamil- 
tonian describes an isolated mesoscopic grain whose single- 
particle dynamics are chaotic (or weakly diffusive in a dis- 
ordered grain) in the limit where the Thouless conductance 
g T — >. oo. It can be derived from general symmetry considera- 
tions^H 

Although the form of the universal Hamiltonian in Eq. <£TJ is 
based on the chaotic (or diffusive) nature of the single-particle 
states, we do not study here the mesoscopic fluctuations, but 
assume a generic equidistant single-particle spectrum (i.e., a 
picket-fence spectrum) as our benchmark model. We consider 
a half-filled band of 2N a + 1 doubly degenerate levels. The 
even grain contains N = 2N electrons while the odd grain 
contains N = 2N + 1 electrons. The single-particle energies 
are given by = kd with k = —N , . . . ,N . All energy scales 
in this work are measured in units of the single-particle mean 
level spacing 8, and for simplicity we take 8=1. 

The reduced BCS model with an attractive pairing force 
is characterized by two regimes: the fluctuation-dominated 
regime or perturbative regime A/8 <C 1 (A is the zero- 
temperature BCS gap), and the BCS superconducting regime 
or non-perturbative regime A/8 ~S> 1. Here we study the ther- 
modynamics of the universal Hamiltonian for three differ- 
ent values of A/8 in the crossover between the fluctuation- 
dominated regime and the BCS regime: A/8 = 0.5, 1 and 5. 
The effective pairing strengths G that correspond to these BCS 
gaps depend on the band width and are calculated using the 
appropriate renormalization methodi 12,25 

The thermodynamic properties of the reduced BCS model 
(in the absence of exchange interaction) are universal func- 
tions of T/5 that depend only on A/8, i.e., changing the 
model-space size for a fixed A/8 and renormalizing G leaves 
the thermodynamic quantities invariant. 12 Of course, choosing 
a smaller model space restricts the temperature range in which 
the model is physically meaningful because of truncation ef- 
fects. In this work we calculate thermodynamic properties for 
even (odd) grains with N = 50 (JV = 51) electrons in a half- 
filled band (of width N = 25) around the Fermi energy in the 
presence of both pairing and exchange correlations. As long 
as the number of blocked levels (i.e., singly occupied levels) is 
much smaller than the total number of levels in the band, the 
thermodynamic quantities are still universal function of T/8, 



but now they depend on two parameters: A/8 and 7 4 /8. As the 
bandwidth is truncated, J s remains invariant while the renor- 
malization of G is approximately independent of J s . We have 
tested this numerically; thermodynamic functions for a band 
width of N = 50 were reproduced by considering a grain with 
a band width of N = 25 and an appropriately renormalized 
coupling strength G. As an example, we show in Fig. Q] the 
heat capacity and spin susceptibility for even and odd grains 
with N = 50 and N = 25. The pairing strength is renormal- 
ized to keep the BCS pairing gap fixed at A/8 = 5, while the 
exchange coupling is fixed at J s = 0.68. 



CL 




FIG. 1 : The heat capacity C and the spin susceptibility % [normalized 
to its bulk high temperature limit %P (Pauli susceptibility) for J s = 0] 
as a function of temperature T for an even and an odd grain. Results 
are shown at half filling for band widths of N a = 25 and N a = 50. 
The pairing strength G is renormalized such that the BCS pairing 
gap is kept fixed at A/8 = 5, while the exchange constant is fixed 
at J s = 0.68. The coincidence of the corresponding thermodynamic 
quantities for both band widths demonstrates that the renormalization 
of G is approximately independent of J s . 



Throughout this work we consider values of the spin cou- 
pling constant J s ranging from to 0.91 8. Values for J s rang- 
ing from / s /8 k -0.03 - 0.09 for copper to 7 v /8 w 0.84 - 0.89 
for palladium were reported in Ref. 26 (extracted from both 
experiment and theory). Since all thermal averages we calcu- 
late are canonical (i.e., for a fixed number of electrons), the 
charging energy term EcN 2 is just an overall constant. We 
therefore put Ec = in our calculations without loss of gen- 
erality. 
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III. QUANTUM MONTE CARLO APPROACH 



IV. THERMODYNAMIC PROPERTIES 



To compute thermodynamic properties of a metallic grain, 
we use a quantum Monte Carlo (QMC) method that is based 
on the canonical loop updates of Refs. I27ll28l There, it was 
shown that this method can be used to simulate the reduced 
BCS Hamiltonian [i.e., Eq. ([D with J s = 0] in the canonical 
ensemble at finite temperature. The QMC method starts from 
a perturbative expansion of the partition function at inverse 
temperature (3 



Z^nJo JO JO 



Tr[y(Ti)y(T 2 )---v(x m ) 



(2) 



where V(l) =exp( — zHo)Vexp(xHo)- The Hamiltonian in 
Eq. (O is assumed to consist of two non-commuting parts, ftp 
and V. In case of the reduced BCS Hamiltonian, these were 
chosen to be 



H = 
Hd = 

V = 
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The basic idea of the QMC method is to insert a so-called 
worm operator A in the partition function, obtaining an ex- 
tended partition function Tr(Ae~$ ). By propagating this 
worm operator th rough i maginary time according to the rules 
explained in Refs. l27ll28l one generates configurations that are 
distributed according to the weights occurring in the canoni- 
cal partition function Tr^yfe - ' 3 ) through a Markov process. 
The worm propagation rules are constructed such that the de- 
tailed balance condition is satisfied. In case of the reduced 
BCS Hamiltonian, the worm operator consists of two parts: 
one that enables scattering of S = pairs, and another that en- 
ables the breakup of an S = pair (thus creating two blocked 
levels). 

To study the universal Hamiltonian, it is necessary to in- 
clude the exchange interaction term. In general, terms that 
commute with ftp can be incorporated in the current algorithm 
by adding them to Ho- The exchange term — J S S 2 commutes 
with Hd in Eq. (0J, and only unpaired electrons (that block 
levels) contribute to the total spin S. For a given number b of 
blocked levels, the degeneracy of many-particle levels in the 
total spin S is given by 



d b (S) 



S- 



b 
1- 



(6) 



Since the number b of blocked levels is known at each step of 
the Markov process, one can simply take the spin exchange 
term into account by adding it to Hp and choosing the total 
spin of the configuration with a probability proportional to the 
degeneracy db(S). The non-diagonal part V remains the same 
as for the reduced BCS model [see Eq. ©I and there is no 
change in the canonical loop updates. 



In the following we use the QMC method to study various 
thermodynamic properties of the grain. 



A. Thermal spin distributions 

We first study the spin distribution at fixed temperature. For 
that purpose, we consider the ratio of the spin-projected parti- 
tion function Z$ (at spin S) to the total partition function for a 
fixed number of electrons 



Z s _ TiNse-W 
Z Trjve-P" ' 



(7) 



Zg is normalized such that Y,s (25+ 1)Z$/Z = 1 [i.e., the (2S + 
l)-fold degeneracy in the spin-projection quantum number M 
is not included in Zs] . 

We first discuss the case of a pure exchange interaction 
(G = 0), for which the ratios (|7]i can be expressed in closed 
form in terms of canonical quantities of non-interacting spin- 
less fermions (G = J s = 0) using the method of Ref. [2(J For 
G = 0, we can rewrite (O as 



Z S _ ^^Tr^-^o 

z ~ I s (25+l)eP^+ 1 )Tr w , s e-P A o ' 



(8) 



where Hq is the non-interacting Hamiltonian. The spin- 
projected quantities can be calculated from the corresponding 
M-projected quantities using 

Tr w , 5e -P«° = Tr N/i=s e-&> ~ Tv NM=s+l e-^" . (9) 

The traces on the r.h.s. of Eq. (O can be evaluated in terms of 
two particle-number projections that correspond to the num- 
ber of spin-up and the number of spin-down electrons. This 
leads tc3 



Tr w , s e 



P^O _ g-P^V/2+Sg _ P^iV/2-S _ e _ P^V/2+S+l e _ P f lv/2-S-l 

(10) 

where F q is the canonical free energy of q non-interacting 
spinless fermions with a single-particle spectrum e^. 

In the presence of both exchange and pairing correlations, 
we evaluate the ratio Zs/Z using the QMC method outlined in 
Sec. |IIl] The configurations generated in the Markov process 
are distributed according to the weights appearing in the par- 
tition function Tr(e _ P") . Since the degeneracy in S is known 
for each configuration [see Eq. ©], the ratio Zs/Z can be 
evaluated directly through 



Z s _/d b (S)\ 
Z \ 2 b Imc 



(11) 



where (...)mc denotes averaging over all the configurations 
generated by the Monte Carlo method, db(S) is the degeneracy 
defined in Eq. (O, and b is the number of blocked levels in the 
configuration. 
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FIG. 2: Thermal ratios fs of the spin projected partition function to 
the full partition function as a function of temperature for an even 
number of electrons (N = 50). The left column shows the fractions 
fs in absence of the pairing interaction (A/5 = 0) for different spin 
couplings J s (shown in units of 8). The middle (right) column corre- 
sponds to to a gap of A/8 = 0.5 (A/8 = 5). The different spin values 
S are indicated by different symbols: o (5 = 0), □ (5 = 1) and A 
(S = 2). 



A/8 = 0.0 A/5 = 0.5 
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FIG. 3: AsinFig.|2]butfor a grain with an odd number of electrons 
(N = 51). The half-integer spin values S are now indicated by the 
following symbols: o (5 = 1 /2), □ (S = 3/2) and A (S = 5/2). 



We define the thermal fraction fs of a given spin 5 by 

/s = (25+l)J. (12) 

Figure[2]shows f s as a function of temperature for a grain with 
an even number of electrons. Results for a few lowest spin 
values are shown for various values of the exchange coupling 
J s (measured in units of 8) and pairing gap A/8. 

The left column of Fig. [2] corresponds to electrons inter- 
acting only through spin exchange (G = and thus A/8 = 
0). In the absence of both pairing and exchange interactions 
(G = J s = 0), the ground state for an even number of electrons 
is found by filling the lowest single-particle energy levels by 



spin up/spin down electrons resulting in an 5 = ground state. 
At low temperatures (T < 0.758), the 5 = states give the 
largest contribution to the partition function. For higher tem- 
peratures, the contribution of higher spin states increases, and 
in the temperature region 0.758 < T < 3.5 8 the largest con- 
tribution arises from the 5=1 states. 

The exchange interaction shifts down in energy states with 
5^0; thus less thermal energy is required to excite these 
states and the fractions fs of non-zero spin values increase 
with J s . For J s = 0.48, the 5 = states dominate only below 
T = 0.28. When J s > 0.58, the ground state acquires a finite 
non-zero spin, and has spin 5 = 3 for J s = 0.845 8. 

Table U lists the values of J s at which the ground-state spin 
changes to a higher value (denoted by S) for A/8 = 0, 0.5 and 
1. For a strong pairing interaction (A/8 = 5), the system re- 
mains fully paired up to an exchange coupling of J s /8 ~ 1, at 
which it makes a transition to a fully polarized state. These 
values of J s were obtained using Richardson's solution to the 
pairing Hamiltonian via the method of Ref. The ground- 
state spin diagram in the presence of pairing correlations and 
ferromagnetism at zero temperature was discussed in Ref.l23i 

In the presence of pairing correlations, low spin states are 
favored because the scattering of spin zero pairs lowers the 
free energy. The middle column of Fig.|2]shows the fraction fs 
for a weak pairing force (A/8 = 0.5). Comparing the fractions 
fs at J s = and J s = 0.48 with their corresponding values in 
the absence of pairing, we observe that the pairing interaction 
makes the 5 = channel more dominant at low temperatures. 
At higher temperatures pairing correlations are destroyed by 
thermal excitations, and there is almost no difference between 
the A/8 = and A/8 = 0.5 cases. The results of a strong pair- 
ing force (A/8 = 5) are shown in the right column of Fig. [2] 
5 = states remain dominant up to higher temperatures and 
the spin fractions are less affected by the exchange interaction. 

As we increase the pairing strength at fixed J s , the 5 = 
channel becomes more dominant at low temperatures, and 
higher values of J s are required to make the transition to a 
higher spin ground state. This in turn affects the finite tem- 
perature behavior of the grain. At a fixed pairing gap A/8 and 
for increasing J s , the crossing point where the 5=1 channel 
becomes dominant shifts to lower temperatures. 

Figure [3] shows results analogous to Fig. [2] but for a grain 
with an odd number of electrons and thus half-integer spin. In 
the absence of pairing, the odd grain has an 5 = 1/2 ground 
state for J s = and acquires higher spin for a sufficiently 
strong exchange interaction. As compared with the even case, 
higher values of J s are required to make the respective transi- 
tions to higher spin states (see Table I). 



B. The number of 5 = electron pairs 

Since our QMC method works directly in the space where 
spin is a good quantum number, we can evaluate the number 
of 5 = pairs for each sampled configuration. Hence we can 
calculate the average fraction of 5 = pairs 
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TABLE I: The J s values (in units of 8) at which the ground state of 
an even grain (left panel) and an odd grain (right panel) acquires a 
higher spin value 5, as J s is increased at fixed A/8. Three values 
of A/8 are considered (0, 0.5 and 1). For A/8 = 1 and N even, the 
ground-state spin makes a transition from S = to S = 4 at J s = 
0.9295 8. 
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where /V the total number of electrons, b the number of 
blocked (i.e., singly occupied) levels and p the parity of the 
grain, i.e., p = (p = 1) for /V even (odd). The normaliza- 
tion in Eq. ( TT3l ) is chosen to give fp — 1 at zero temperature 
for both the even and odd grain in the absence of spin ex- 
change. Figure [4] shows fp as a function of temperature for 
even (open circles) and odd (solid triangles) grains with BCS 
gaps of A/8 = 0.5, 1 and 5. In general, the number of pairs 
decreases with temperature, reflecting the weakening of pair- 
ing effects with increased thermal energy. The temperature 
at which pairs start to break up depends strongly on A/8. In 
the case of strong pairing with A/8 = 5 this temperature is 
roughly ~ 1 .5 8 (for J s — 0) and reduces to a value of ~ 0.48 
in the weak pairing case with A/8 = 0.5. 

Even as pairs start to break, their number decreases only 
slowly with increasing temperature and for T ~ 3.5 8 most 
of the electrons are still paired to S — (e.g., about 84% for 
A/8 = 1 and J s = 0.6). Once pairs start to break up, the frac- 
tion of pairs fp is always larger for the odd grain. This is be- 
cause the extra electron blocks a level, deferring the transition 
to higher temperatures. Consequently more thermal energy is 
required to break up the same number of pairs in the odd grain 
as compared with the even grain. 

The overall effect of the exchange interaction is to reduce 
the average number of pairs. Exchange also reduces the 
threshold temperature at which pairs start to break. For a pair- 
ing gap of A/8 = 1 this temperature is about ~ 0.58 for J s = 0, 
and it reduces to ~ 0.28 for J s = 0.8 8. Indeed, in an even 
grain the exchange interaction decreases the gap between the 
5 = ground state and the first 5^0 excited state, thereby 
reducing the thermal energy required to break up a pair. For 
A/8 = 0.5 and J s = 0.845 8, the ground state of the even grain 
is S = 1 so there are two blocked levels at T = 0. 

Pairing correlation effects can be more clearly observed in 
the increase of the number of S = pairs as we turn on the 
pairing interaction at a fixed exchange interaction, i.e., (n p ) — 
( n p)c=o at fixed J s . This pair number excess is simply related 
to the deficiency of the average number of blocked levels 

(n P )-(n p ) G= o= 1 -{(b) G=0 -(b)) . (14) 

Figure|5]shows the excess number of S — pairs for both even 
and odd grains. A clear odd-even effect is observed in this 



FIG. 4: The fraction fp of 5 = pairs as a function of temperature for 
pairing gaps of A/8 = 0.5 (left column), A/8 = 1 (middle column) 
and A/8 = 5 (right column), and for different values of the exchange 
coupling J s (shown in units of 8). Results for the even grain are 
shown by circles (o), while results for the odd grains are denoted by 
solid triangles (A). 



quantity. For a stronger pairing interaction (larger A/8), the 
odd-even effect survives up to higher temperatures. For not 
too large J s , the excess number is zero at T = because the 
ground state has the minimal spin and therefore the largest 
possible number of pairs even in the absence of pairing inter- 
action. For the largest values of J s (third row in Fig. |5]l, the 
excess number of pairs is non-zero at T = since the grain is 
polarized. 



A/8 = 0.5 A/8 = 1 A/8 = 5 



CD 



J s = 0.0 


J s = 0.0 


J s = 0.0 


J s = 0.4 


J s = 0.4 


fx 

'■"f iii 


J s = 0.845 


J s = 0.8 





0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 2.5 



T/S 

FIG. 5: The excess number of 5 = pairs [see Eq. dl4H as a function 
of temperature. We show results for both even (o) and odd (A) grains 
with A/8 = 0.5 (left column), A/8 = 1 (middle column) and A/8 = 5 
(right column). 
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C. The canonical pair gap 



D. The heat capacity 



The canonical pair gap A can , defined by 



A c 2 an (r,G,/ ( ) =G((P f P) G ,j s - (P t P} G =oj 



Another interesting thermodynamic observable is the heat 
capacity of the grain 



(15) 



measures the pairing correlation energy, namely the increase 
of pairing energy when the pairing interaction is turned on in 
the presence of a fixed exchange interaction. For J s = and in 
the thermodynamic limit, the canonical pair gap A can becomes 
the familiar BCS gap A. For a finite system, the BCS gap A 
is recovered from A can by applying the mean-field approxima- 
tion and taking the grand-canonical averages in Eq. (TT3T >. 

We first discuss the behavior of the canonical pair gap for 
a weak pairing interaction A/8 = 0.5 (left column of Fig. [6]). 
We observe that the pairing correlation energy decreases with 
increasing temperature. The behavior of A can versus temper- 
ature is completely smooth because of the finite size of the 
grain. The exchange interaction quenches the pairing corre- 
lation energy further since this interaction tends to break up 
pairs. 

At low temperatures and J s not too large, an odd-even dif- 
ference is visible in the canonical pair gap. This odd-even ef- 
fect is a unique signature of pairing correlations and is reduced 
by the exchange interaction. For J s = 0.845 8 the ground-state 
spin of the even (odd) grain is S = 1 (5 = 3/2) and the odd- 
even effect in A can is completely destroyed. 

For larger pairing strengths, the exchange interaction does 
not affect much the canonical pair gap, as can be seen from the 
middle and right columns of Fig. [6] corresponding to A/S = 1 
and A/8 = 5, respectively. 



A/8 = 0.5 



A/5= 1 



A/5 = 5 




FIG. 6: The canonical pair gap A can defined in Eq. U5\ as a function 
of temperature. We show results for both even (o) and odd (A) grains 
with A/8 = 0.5 (left column), A/8 = 1 (middle column) and A/8 = 5 
(right column). Visible effects of the exchange interaction are limited 
to the weak pairing case (A/8 = 0.5). 



c 



d(H) 
dT 



(16) 



with (...) denoting thermal averaging. Figure [7] shows the 
heat capacity in grains with BCS gaps of A/8 = 0.5, 1 and 5. 

We first discuss the smaller grains with a BCS gap of 
A/8 = 0.5 (left column of Fig. |7). Previous studies have 
shown that in the absence of exchange interaction (J s = 0) the 
even-grain heat capacity exceeds the odd-grain heat capacity 
in a temperature range 0.48 < T < 1 .3 8 . 11 ' 12 In this temper- 
ature range, S = pairs start to break up and pairing corre- 
lations are quenched, as can be seen from Fig. [6] The bump 
in the heat capacity of the even grain reflects a signature of 
the pairing transition of the finite-size grain. No such effect 
is observed in the odd case because of the blocking effect of 
the unpaired electron. This odd-even effect in the heat capac- 
ity is a unique signature of pairing correlations in a finite-size 
system. 

When the exchange interaction is turned on, the odd-even 
effect in the heat capacity disappears gradually (see left col- 
umn of Fig. |7j. Since the exchange interaction brings down 
in energy high spin states while leaving the 5 = states un- 
changed, it increases the number of unpaired electrons at finite 
temperature (even when the even ground state still has 5 = 0). 
These unpaired electrons block levels (in the same way as the 
single electron blocks a level in the odd grain) and suppress 
the bump in the heat capacity. For J s = 0.8458, the ground 
state has S = 1 (S = 3/2) in the even (odd) grain. For this 
value of J s , the odd-even effect has completely disappeared. 
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FIG. 7: The heat capacity for even (o) and odd (A) grains as a func- 
tion of temperature. Shown are results for grains with A/8 = 0.5 (left 
column), A/8 = 1 (middle column) and A/8 = 5 (right column). 



A similar behavior is found for grains in the crossover re- 
gion A/8 = 1 (middle column of Fig.[7]l, where the signature 
of the pairing transition is destroyed by the exchange inter- 
action. Compared to the smaller grains (with A/8 = 0.5), 
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the odd-even effect is larger, and a larger critical exchange 
strength is required to destroy it. 

The right column of Fig. [7] shows the heat capacity of a 
grain in the BCS regime (A/8 = 5). The qualitative differ- 
ence with the behavior of the heat capacity in the fluctuation- 
dominated regime (A/8 = 0.5) is striking: the signature of the 
pairing transition is much stronger and it cannot be destroyed 
even in the presence of a strong exchange interaction. We 
can understand this effect by comparing the fraction of 5 = 
pairs in both cases (see left and right column of Fig. |4j». When 
the exchange interaction is increased, it is clear that 5^0 
states are pushed down in energy. At some critical value of 
the exchange strength, the ground state eventually acquires a 
finite spin. Before this happens, the gap between the 5 = 
ground state and the first 5^0 excited state decreases with J s . 
However, due to strong pairing correlations, this gap is much 
larger in the BCS regime. Thus, in the fluctuation-dominated 
regime, little thermal energy is needed to excite the system 
to 5 ^ states, whereas in the BCS limit 5 = states dom- 
inate up to considerably higher temperatures. This effect is 
also reflected in the number of 5 = pairs (see Fig. |4). For 
A/8 = 5, the excitation gap is still large enough and a clear 
finite temperature transition occurs at a considerably higher 
temperature (~ 1.5 8) from a 5 = state to a state with bro- 
ken pairs due to thermal excitations. For this A/8 = 5 case, 
the (even) system makes a sudden transition at J s = 1 .0029 8 
from a 5 = ground state to a state where all the electrons 
in the model space are unpaired. This is known as the Stoner 
instability^ 



E. The spin susceptibility 

The spin susceptibility is a measure of the grain's response 
to an external magnetic field. Here we discuss the spin sus- 
ceptibility in the zero-field limit, defined by 



A/5 = 0.5 



A/8= 1 



A/5 = 5 



X(T) 



d 2 f(T,h) 



dh 2 



a4 



h=0 



«M 2 ) - (M) 7 



(17) 




FIG. 8: The spin susceptibility, normalized to its J s = bulk high 
temperature limit %p = 2/^/5, as function of temperature for an even 
(o) and odd (A) grain with A/5 = 0.5 (left column), A/5 = 1 (middle 
column) and A/5 = 5 (right column). The ground-state spin is S = 
or S = 1/2 for all shown values of the spin coupling J„ (measured in 
units of 5). The dotted lines are the even and odd spin susceptibility 
for electrons interacting only through an exchange interaction with 
the indicated strength J s (in units of 5). 



step of the Markov process. To quantify the effects of pairing 
correlations, we compare our results with the limiting case of 
spin exchange correlations but no pairing interaction. In this 
limiting case, the spin susceptibility can be calculated directly 
using spin projection methods. 20 We find 



X(T) 



I s (25+l)e-^%P- / ^+ 1 ) 



(18) 



where 



and 



where J (T,h) = -Tln[TTe-^-^" M )] is the free energy of 
the grain in the presence of an external Zeeman field h and 
g is the spin g-factor. The operator M is the "magnetization" 

defined as M = L,', C7cJ a c ( . c . For h = 0, we have (M) = be- The quantity F q is the canonical free energy of q non- 



-4(5+ l) 2 e -P( f 'V/2+5'+l+ f lv/2-S-l) ] (19) 
,~$ F N,S — e —$( F N/2+S+*N/2-S) _ e _ P(^V/2+5+l+^iV/2-S-l) _ (20) 



cause of spherical symmetry. Within the reduced BCS model, 
it was found that pairing correlations affect the temperature 
dependence of the spin susceptibility of a grain i 9 ' 11,12 In par- 
ticular, for an odd number of electrons, the spin susceptibility 
shows a re-entrant behavior as a function of T for any value 
of the ratio A/8. This behavior persists in ultra-small grains, 
in which the level spacing is larger than the BCS gap. Since 
this re-entrant behavior is absent in normal metallic grains, it 
was suggested by Di Lorenzo et al. 9 that this behavior could 
be used as a unique signature of pairing correlations in small 
grains. Here we study how the exchange interaction affects 
this re-entrant behavior. 

It is straightforward to evaluate the spin susceptibility in 
the QMC approach since the value of M is known at each 



interacting spinless fermions in 2N + 1 single-particle levels, 
which can be evaluated using a particle number projection for- 
mula that involves 2N + 1 quadrature points*^ 

The left column of Fig. [8] shows the even and odd spin 
susceptibility for weak pairing (A/8 = 0.5) and for spin ex- 
change couplings J s < 0.68. For these values of J s , the even 
(odd) ground state has 5 = (5 = 1/2) (see Table IT]). The 
susceptibility is measured in units of the Pauli susceptibility 
%p = 2/j|/8 (the high-temperature value of % at J s = 0). In 
particular, the top left panel of Fig. [8] shows the spin suscep- 
tibility in the absence of exchange interaction (J s = 0). At 
low temperatures, the spin susceptibility is exponentially sup- 
pressed for the even grain, but exhibits the familiar re-entrant 
effect for the odd grain. This re-entrant behavior is seen for 
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all cases with A/8 = 0.5 and J s < 0.68, including the case 
J s = 0.6 8 for which the signature of pairing correlations is no 
longer visible in the heat capacity. 

This re-entrant behavior originates in the paramagnetic con- 
tribution of the spin of the unpaired electron. This contribu- 
tion is given by %(T)/%p = 5/2T (not shown in the figure), 
and coincides with the odd-grain QMC results at sufficiently 
low temperatures (T < 0.48). At higher temperatures, the 
QMC results deviate from this simple behavior since several 
unpaired electrons contribute to the odd-grain susceptibility. 
These deviations are correlated with the breakup of S = pairs 
(see Fig. @). The stronger the pairing strength, the higher the 
temperature at which the spin susceptibility deviates from the 
b/2T behavior, as more thermal energy is required to break 
up pairs. 

For comparison, we also show in Fig. [8] the even and odd 
spin susceptibilities when the electrons interact only through 
the exchange channel (dotted lines). In general, we observe 
that exchange correlations enhance the spin susceptibility. For 
J s = 0.48 and A/8 = 0.5, we observe a peak in the even spin 
susceptibility (in the absence of pairing) around T w 0. 1 8. For 
this value of J s , the excited triplet state (5=1) lies close to the 
ground-state singlet (S = 0) and the system could be easily 
polarized at low temperatures (the S = — >• S = 1 ground- 
state spin transition occurs at J s = 0.58). At temperatures T < 
0. 1 8, the even spin susceptibility is exponentially suppressed, 
while at slightly higher temperatures the susceptibility tends 
to follow the odd spin susceptibility. At J s = 0.68, both the 
even and odd spin susceptibility diverge at T = 0, since the 
ground state has already acquired a finite spin (for A/8 = 0). 

We also observe from Fig. [8] that at high temperatures the 
spin susceptibility in the presence of pairing correlations ap- 
proaches its value in the absence of pairing. This behavior 
is expected since pairing correlations are suppressed at high 
temperatures. 

The middle column of Fig.[8]shows the spin susceptibilities 
for even and odd grains with BCS gap of A/8 = 1. For the 
largest exchange value shown (J s = 0.8 8) the ground state is 
still S — (or S — 1 /2). At low temperatures we observe (for 
J s = 0.88) a clear minimum in the odd spin susceptibility (a 
signature of pairing correlations), and the 8/T behavior of a 
single unpaired spin is observed now only for T < 0.1 8. For 
T > 8, we note the decrease of the spin susceptibility with 
temperature. 

The case of strong pairing with A/8 = 5 is shown in the 
right column of Fig. [8] S = pairs start to break up only at 
higher temperatures and the behavior of a single spin suscep- 
tibility of ~ 8/T for the odd grain remains valid up to T « 8. 
At high temperatures, the spin susceptibility increases with 
the exchange coupling. At this large pairing strength, the first 
spin jump occurs at J s = 1.00298, and it immediately polar- 
izes the entire system. 

Figure [9] shows the spin susceptibility for A/8 = 0.5 and 
spin exchange values of J s > 0.8 8. At J s = 0.88, the system 
is close to its first ground-state spin jump (which occurs at 
J s = 0.83798 for an even grain or at J s — 0.83208 for an odd 
grain), and the system is easily polarized. The exponential 
suppression of the even susceptibility can only be observed at 
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FIG. 9: The spin susceptibility as function of temperature for an even 
(o) and odd (A) grain with A/8 = 0.5. For spin coupling values J s 
of 0.88, 0.8458, 0.8858 and 0.918, the even (odd) ground state has 
spins 5 = (5 = 1/2), S = 1 (5 = 3/2), 5 = 2 (S = 5/2) and S = 3 
(5 = 7/2), respectively. 



very low temperatures (T < 8/40). The susceptibility peaks 
at T « 0.058. At higher temperature, there is a large num- 
ber of broken pairs and the even spin susceptibility coalesces 
with the odd susceptibility. The odd spin susceptibility is a 
monotonic function and no re-entrant behavior is observed. 

Once the ground-state spin transition has occurred (J s > 
0.8458), both the even and odd susceptibilities diverge at 
T = with the even curve lying slightly below the odd curve. 
For J s = 0.845 8, the odd-grain ground state has S = 3/2, and 
we expect a low temperature behavior of %{T) /%p = 58/2T 
(dashed line). For J s = 0.8858 (7,. = 0.91 8), the ground state 
has S = 5/2 (S = 7/2), leading to a low temperature behavior 
of %(T)/Xp = 358/67 (x(T)/xp = 218/27) of the odd spin 
susceptibility. 

We conclude that once the exchange strength gets close to 
its value where the first ground-state spin jump occurs, the re- 
entrant behavior in the odd spin susceptibility disappears. We 
emphasize, however, that pairing correlations still exist since 
the canonical pair gap does not vanish. 



V. CONCLUSION 

We have used a quantum Monte Carlo method to calcu- 
late the thermodynamic properties of a small superconducting 
metallic grain that is described by the universal Hamiltonian. 
These thermodynamic properties have been studied as a func- 
tion of the BCS gap A/8 and the exchange interaction strength 
J s /8 (measured in units of the mean-level spacing). The spin 
exchange interaction competes with the BCS-like pairing in- 
teraction, and, in general, we find that number-parity signa- 
tures of pairing correlations are suppressed in the presence 
of a finite exchange interaction. We also find qualitative dif- 
ferences between the superconducting BCS regime and the 
fluctuation-dominated regime of pairing correlations. 
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